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OBJECTIVES 

The  objectives  of  this  project  were  to  develop  new  hybrid  analytical/computational 
methods  that  are  capable  of  simulating  the  rare  events  that  are  the  determining  factors  of 
the  performance  of  lightwave  systems  and  devices.  These  methods  use:  1)  Analytical 
techniques,  such  as  perturbation  and  asymptotic  methods,  to  guide  numerical  simulations 
using  importance  sampling;  2)  Adaptive  numerical  methods,  such  as  the  multicanonical 
Monte  Carlo  and  cross-entropy  methods,  to  perform  the  simulation  of  rare  events  when 
guiding  analytical  models  are  not  available.  The  above  methods  can  be  used  to  evaluate 
the  performance  of  specific  optical  systems  and  devices,  including  ultra-high-precision 
optical  clocks  based  upon  mode-locked  fiber  lasers;  optical  clocks  and  other  devices 
based  upon  hybrid  opto-electronic  oscillators.  In  each  case,  the  goal  is  to  use  the 
methods  to  develop  models  that  can  accurately  predict  the  performance  of  these  devices, 
as  well  as  determine  the  failure  modes  that  are  the  limiting  factors  in  their  performance. 


STATUS  OF  EFFORT 

Importance  sampling  has  been  applied  to  compute  the  performance  of  a  differential- 
phase-shift  keyed  (DPSK)  system.  The  results  of  basic  method  were  published  in 
Photonics  Technology  Letters,  and  more  detailed  publications  explaining  the  methods 
have  been  published  in  the  Society  for  Industrial  and  Applied  Mathematics  (SIAM)  Journal 
on  Applied  Mathematics,  and  submitted  to  the  SIAM  Journal  on  Applied  Dynamical 
Systems.  The  method  has  been  optimized,  and  a  numerical  implementation  ported  to  run 
on  a  64-processor  Beowulf  cluster.  This  implementation  was  used  to  determine  the 
optimal  placement  of  an  in-line  phase-conjugator  (the  conjugator  improves  system 
performance  by  reducing  phase  noise).  This  work  was  presented  at  the  Conference  on 
Lasers  and  Electrooptics  in  Baltimore,  Maryland.  We  next  applied  the  importance- 
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sampling  method  to  the  simulation  of  a  mode-locked  laser.  The  method  allows  amplitude 
drop-outs  and  position  shifts  to  be  efficiently  simulated  in  these  systems.  This  work  was 
reported  at  the  2007  Conference  on  Lasers  and  Electro-Optics,  held  in  Baltimore, 
Maryland.  Finally,  a  new  iterative  numerical  method  employing  the  singular  value 
decomposition  and  based  on  the  cross-entropy  method  has  been  developed. 


ACCOMPLISHMENTS/NEW  FINDINGS 

We  have  developed  methods  based  upon  soliton  perturbation  theory  and  importance 
sampling  to  simulate  rare  events  in  lightwave  systems,  including  mode-locked  laser 
systems.  A  key  step  to  using  the  methods  based  upon  soliton  perturbation  theory  is  to 
use  an  approximate  version  of  the  system  dynamics  to  determine  the  locations  in  the 
large-dimensional  state  space  that  most  contribute  to  the  desired  rare  events  (e.g., 
errors).  In  this  method,  calculus  of  variations  applied  to  the  approximate  system  allows 
the  most  significant  rare  events  to  be  located,  and  then  fully-detailed  importance-sampled 
Monte-Carlo  simulations  in  the  vicinity  of  these  locations  properly  determines  the 
probabilities  of  these  rare  events  and  corrects  for  any  errors  made  by  the  approximations 
in  determining  the  system  dynamics.  In  the  case  of  the  simulations  of  mode-locked 
lasers,  the  analysis  shows  that  the  system  dynamics  are  related  to  the  classical  exit-time 
problem  of  stochastic  differential  equations.  For  example,  if  a  periodic  signal  is  added  to 
control  the  pulse  position,  pulses  tend  to  locate  themselves  at  the  bottom  of  an  effective 
potential  well  created  by  the  periodic  signal,  and  noise  causes  the  pulses  to,  once  in  a 
great  while,  to  climb  up  the  potential  well  and  escape  into  a  neighboring  potential  well.  If 
the  noise  is  weak  the  probability  of  these  rare  events  is  very  small. 

More  recently  we  have  developed  a  new  fully-numerical  method  based  upon  the  singular 
value  decomposition  and  the  cross-entropy  method.  The  singular  value  decomposition  is 
applied  to  the  equations  linearized  about  a  specific  nonlinear  solution  and  allows  one  to 
determine  numerically  the  specific  perturbations  that  produce  the  largest  change  in  some 
desired  output  quantity  or  performance  measure,  and  thus  replaces  the  step  of 
determining  the  approximate  dynamics  using  soliton  perturbation  theory.  Because  it  is 
fully  numerical,  non-soliton  pulses  can  be  examined.  Next,  the  most  probable 
perturbations  are  computed  using  the  cross-entropy  method.  This  method  uses  the 
Kullback-Liebler  probability  distance  measure  to  compute  the  deviation  between  the 
current  Monte-Carlo  sampling  distribution  and  the  optimal  importance-sampling  biasing 
distribution.  By  iterating  to  decrease  the  difference  between  the  current  and  the  optimal 
distribution  (which  can’t  be  computed  directly  itself),  one  can  reach  the  most  probable 
locations  in  state  space  that  produce  errors.  Because  the  techniques  used  in  this  method 
are  very  general,  we  anticipate  that  it  can  be  used  to  explore  nonlinear  systems  for  which 
rare  events  are  significant  other  than  lightwave  systems. 
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